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Abstract. I describe the recent measurement of the CP violating parameter sin 2(|)i = (0.99±0.14± 
0.05) from the Belle experiment at KEK. 


INTRODUCTION 


In 1973, Kobayashi and Maskawa (KM) first proposed a model where CP violation is 
ineorporated as an irredueible eomplex phase in the weak-interaetion quark mixing ma¬ 
trix [0- The idea was remarkable and daring beeause it required the existence of six 
quarks at a time when only the u, d and s quarks were known to exist, The subsequent 
discoveries of the c, b and t quarks, and the compatibility of the model with the CP viola¬ 
tion observed in the neutral kaon system led to the incorporation of the KM mechanism 
into the Standard Model, even though it had not been conclusively tested experimentally. 

In 1981, Sanda, Bigi and Carter pointed out that a consequence of the KM model 
was that large CP violating asymmetries could occur in certain decay modes of the 
B mesons. [0] These asymmetries may occur when a neutral B meson decays to a CP 
eigenstate. In this case the amplitude for the direct decay interferes with that for the 
process where the B meson first mixes to a B meson that then decays to the CP eigenstate. 

The unitarity of the CKM matrix implies that the existence of three measurable 
phases. In the “Nihongo” convention, these are denoted 

). ( 1 ) 

while at SLAC these angles are usually referred to as (3, a and y, respectively. 

A non-zero value of ^cp results in the time dependent asymmetry 


RiB^->fcp)-R{B^^fcp) 

^ R{B^^fcp)+R{B^^fcp) 


^f sin 2^cp ■ sin(Am ■ (f 2 ± f i)), 


( 2 ) 


where is the CP eigenvalue (±1), Am denotes the mass difference between the two 
mass eigenstates and ti and t 2 are the proper time for the tagged-B and CP eigenstate 
decays, respectively. The -I- sign corresponds to the case where the B® and B® are in an 
even L orbital angular momentum state, the — sign obtains for odd L states such as the 
T(45). A determination of Ay thus provides a measurement of sin2(|)cp. 

We note that due to the restrictions of quantum mechanics, time integrated asymme¬ 
tries at the T(45) resonance (which corresponds to the — sign in the above equation) 
are identically zero. Therefore, one must make time dependent measurements. Since the 






pairs of B mesons are produeed nearly at rest in the usual arrangement at threshold, the 
T(45) eenter of mass frame must be boosted. This is aeeomplished by the use of beams 
with asymmetric energies. For example at KEK-B, (3y ~ 0.43, and as a result the typical 
B meson decay length is dilated from 20/jm to about 200/rm, which is measurable with 
double-sided silicon strip vertex detectors close to the interaction point. 

The measurement therefore requires: 

• a large sample of reconstructed B (cc)^° eigenstate decays; 

• a determination of the flavor of the accompanying B (“tagging”); 

• a measurement of Az, the vertex separation between the CP eigenstate and flavor 
tag decays; and 

• a fit to the flavor-tagged vertex distribution to extract sin2(|)i. 

The KEKB high luminosity double storage ring facility was commissioned with re¬ 
markable speed starting in late 1998. This accelerator facility allows us to satisfy the first 
requirement. By the summer of 2001, 29.1 fb^^ was integrated on the T(45). This data 
sample was used for the sin2(|)i measurement and corresponds to 31.3 million BB pairs. 
KEKB uses a ±11 mrad crossing angle to separate the incoming and outgoing beams 
and minimize parasitic collisions. So far no special limitations associated with the cross¬ 
ing angle have been observed (e.g. synchrobetatron oscillations). KEKB now routinely 
achieves peak instantaneous luminosities above 5 x jcm^ jsec with acceptable ex¬ 

perimental backgrounds and trigger rates in the Belle detector. The beam currents are 
still far below the design values and therefore there is room for further improvements in 
luminosity. Much larger data samples are expected in the near future. 

The Belle detector has good lepton identification and high efficiency for both charged 
and neutral particles. It allows the CP eigenstate decays to be efficiently reconstructed. 
Belle is a large-solid-angle magnetic spectrometer that consists of a three-layer silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), a mosaic of aerogel 
threshold Cerenkov counters (ACC), time-of-flight scintillation counters (TOE), and 
an array of CsI(Tl) crystals (ECE) located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux-return located outside of the coil is 
instrumented to identify muons and (KEM). The detector is described in detail 
elsewhere [|^]. Examples of its performance are given in the following section. 

The Belle experiment starting physics data taking in 1999. In the summer of 2001, 
Belle (together with BABAR) announced the observation of the first statistically signif¬ 
icant signals for CP violation outside of the kaon system. In this report I will describe 
some of the details of the measurement. 


SIN 2^1 MEASUREMENT 

We reconstruct 5® decays to the following CP eigenstates [|5|]: J/'^Ks, \[/(25)K5, %c\Ks, 
ricA ^5 for — 1 and J for = ±1. We also use ^ decays where 

K*Q KsTl^. Here the final state is a mixture of even and odd CP, depending on the 
relative orbital angular momentum of the 7/\|/ and K*^. The CP content is determined 
from a fit to the full angular distribution of all J/xfK* decay modes other than > 
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FIGURE 1. The beam-energy constrained mass distribution for all decay modes combined other than 
J/\\fKL. The shaded area is the estimated background. The signal region is the range 5.27 — 5.29 GeV/c^. 


Ks'kP. We find that the final state is primarily = -1-1; the — —I fraetion is 0.19 ± 
0.04(stat) ±0.04(syst) [0. 

J/\\f and tj/(25) mesons are reeonstrueted via their deeays to {£ = /j,e). The 
is also reeonstrueted via and the x^i via J/'^y. The Tjc is deteeted in 

the K^K 'nP and KsK modes. For the J/'^Ks mode, we use Ks and 

deeays; for other modes we only use Ks 7t+7t^. 

The 7/\[f and tj/(25) —> eandidates are reeonstrueted from oppositely eharged 

traek pairs where at least one traek is positively identified as a muon by the KLM system 
and the other is either positively identified as a muon or has an ECL energy deposit 
eonsistent with that of a minimum ionizing partiele. For deeays, we use oppositely 
eharged traek pairs where at least one traek is a well identified eleetron and the other 
traek satisfies minimal dE/dx and E/p requirements. For dieleetrons, we eorreet for 
final state radiation or bremsstrahlung and reeover additional \|/ eandidates by ineluding 
the four-momentum of every photon deteeted within 0.05 radians of the original or 
direetion in the invariant mass ealeulation. Candidate Ks —> deeays are 

oppositely eharged traek pairs that have an invariant mass between 482 and 514 MeV/c^, 
whieh eorresponds to ±3a around the Ks mass peak. 

To reeonstruet XciKs deeays, we seleet %ci ^ deeays, rejeeting y’s that are 
eonsistent with yy deeays, and impose the requirement 385 < < 

430.5 MeV/c^. For Tjc deeays, we distinguish kaons from pions using a eombination 
of CDC dE/dx measurements and information from the TOF and ACC systems. Candi¬ 
date Tjc — K^K^n^ (KsK^n^) deeays are seleeted with a KKn mass requirement that 
takes into aeeount the natural width of the ri^. For J/\\fK*^{Ks'nP) deeays, we use Ks'kP 
eombinations that have an invariant mass within 75 MeV/c^ of the nominal K* mass. 
We reduee baekground from low-momentum 7t®’s by requiring eos 0a:* < 0.8, where 0a:* 
is the angle between the Ks momentum veetor and the K*^ flight direetion ealeulated in 
the K*^^ rest frame. 












TABLE 1. The numbers of observed events 
(Aev) and the estimated background (Abkgd) in 
the signal region for each fcp mode. 


Mode Aev iVbkgd 


J)Ks (Tt+Tt^ ) 

457 

11.9 

J)Ks 

76 

9.4 

\(/(25) ^) Ks (Tt+Tt^) 

39 

1.2 


46 

2.1 

Xci{J/^y)Ks{Ti+n-) 

24 

2.4 

■r\c{K+K-K^)Ks{K+K-) 

23 

11.3 

y\c{KsK-n+)Ks{%+%^) 

41 

13.6 


41 

6.7 

Sub-total 

747 

58.6 

jMe+i-)KL 

569 

223 


Reconstructed B meson decays are identified using the beam-constrained mass 

^beam = ^“ Pb ^hc energy difference AE = Eg - Ebeam, where Eteam is 

the cms beam energy, and ps and Eg are the B candidate three-momentum and energy 
calculated in the cms. 

Figure shows the combined Mbc distribution for all channels other than J 
after a mode-dependent requirement on AE. The B meson signal region is defined as 
5.270 < Mbc < 5.290 GeV/c^. Table lists the numbers of observed candidates (A^ev) 
and the background (A^bkgd) determined by extrapolating the rate from the AE vs. Mbc 
sideband region into the signal region. About 65% of the exclusive CP signal events 
(44% of the full CP sample) are reconstructd in the 5® ^ ^Ks, Ks 7t+7t^ mode. 

Candidate 5® ^ J/\\fKi decays are selected by requiring ECL and/or KLM hit pat¬ 
terns that are consistent with the presence of a shower induced by a neutral hadron. Kt 
candidates with ECL information only are treated separately from the candidates 
with KLM hits. The centroid of the shower is required to be within a 45° cone centered 
on the Ki^ direction that is inferred from two-body decay kinematics and the measured 
four-momentum of the 7/t|/. To reduce the background we cut on a likelihood ratio that 
depends on the J/\\f cms momentum, the angle between the Kl and its nearest-neighbor 
charged track, the charged track multiplicity of the event, the extent to which the event is 
consistent with a 5+ ^ J(Kin^) hypothesis, and the polar angle with respect to 
the z direction of the reconstructed 5® meson in the cms. In addition, events that were re¬ 
constructed as RO J/y\fKs, B+ J/\\fK+, or J/\\fK*+{K+%^, 

KsH^) decays are removed. Einally, clusters with positions that match photons from 
reconstructed 7t®’s are also rejected. 

Eigure ^ shows the distribution, calculated with the 5° — J/\\fKi two-body 
decay hypothesis. The histograms are the results of a fit to the signal and background 
distributions. The shapes are derived from the Belle GEANT based Monte Carlo (MC) 
simulation. However, the normalization and peak position of the signal are allowed to 
vary. There are 397 entries in the 0.2 < pg™® < 0.45 GeV/c signal region with KLM 
clusters. There are 172 entries in the range 0.2 < pg™® < 0.40 GeV/c with clusters in 
the ECL only. The fit finds a total of 346 ± 29 J/\\fKi signal events, and a signal purity 






of 61%. Thus, about 33% of the signal in the full CP eigenstate sample is reeonstrueted 
in the 5® MfKi mode. 



FIGURE 2. The distribution for > J/MfKi candidates with the results of the fit. The solid line 
is the signal plus background; the shaded area is background only. The signal region for KLM (ECL-only) 
clusters is 0.2 < < 0.45(0.40) GeV/c. 

To identify the flavor of the aeeompanying B meson, leptons, kaons, A’s, eharged 
slow pions from D* —^ deeays, and energetie pions from two-body B deeay 

(e.g. 5° are used. A likelihood based method, deseribed in detail below, 

is used to eombine information from the different eategories and to take into aeeount 
their eorrelations. The figure of merit for flavor tagging performanee is the effeetive 
effieieney, £gyy, whieh is £(1 — 2w^) summed over all tagging eategories. This method 
gives Zeff = 0.270 ± O.OOSlJg. 

A Monte Carlo simulation is used to determine a table for a eategory-dependent 
variable that indieates whether a partiele originates from a 5® or BP . The values of 
this variable range from — 1 for a reliably identified 5® to +1 for a reliably identified 
bP and depend on the tagging partiele’s eharge, em momentum, polar angle, partiele- 
identifieation probability, as well as other kinematie and event shape quantities. For 
lepton tags, the missing momentum and reeoil momentum are ineluded in the likelihood 
determination. For slow pion tags, the angle between the pion and the thrust axis of 
the non-fcp traeks is used. Charged kaon tags aeeompanied by Ks mesons, whieh have 
additional strange quark eontent and a lower tagging value, are treated as a separate 
tagging eategory. 

The results from the separate, partiele-level eategories are then eombined in a seeond 
stage that takes eorrelations for the ease of multiple partiele-level tags into aeeount. 
This seeond stage determines two event-level parameters, q and r. The first, q, has the 
diserete values ^ = -I- 1 when the tagged B meson is more likely to be a bP and — 1 when 
it is more likely to be a 5®. The parameter r is an event-by-event flavor-tagging dilution 
faetor that ranges from r = 0 for no flavor diserimination to r = 1 for unambiguous flavor 
assignment. The value of r is used only to sort data into six intervals of flavor purity; the 
wrong-tag probabilities that are used in the final CP fit are determined from data. 



TABLE 2. The event fractions (//) and in¬ 
correct flavor assignment probabilities (w/) 
for each r interval. The errors include both sta¬ 
tistical and systematic uncertainties. 


/ 


r 

ft 

Wl 

1 

0.000 

-0.250 

0.405 

U.H-DJ_qqq9 

2 

0.250 

-0.500 

0.149 

u.:)jz_QQ^4 

3 

0.500 

-0.625 

0.081 

u.zH-:)_0 03o 

4 

0.625 

-0.750 

0.099 

0 176+0-022 
/O_0.017 

5 

0.750 

-0.875 

0.123 

0 11 

U.iiU_ooi4 

6 

0.875 

-1.000 

0.140 

0 041 

U.U^i_0,010 



FIGURE 3. Af distributions for the events with q^f = -fl (solid points) and q^f = —1 (open points). 
The results of the global fit (with sin2(|)i = 0.99) are shown as solid and dashed curves, respectively. 


To avoid dependence on the MC, the probabilities of an incorrect flavor assignment, 
w/ (/ = 1,6), are determined directly from the data for six r intervals using exclusively 
reconstructed, self-tagged.6® andJ/\\fK*^{K^'K^) decays. 

The Z^-flavor of the accompanying B meson is assigned according to the flavor-tagging 
algorithm described above. The exclusive decay and tag vertices are reconstructed using 
the same vertexing algorithm that is used in the CP fit. The values of w/ are obtained from 
the amplitudes of the time-dependent mixing oscillations: (A^of — A^sf)/(A^of + 
A^sf) = (1 ~ 2w/) cos(AmrfAf). Here Aqf and A^sf are the numbers of opposite flavor and 
same flavor events. The value of Am^ is fixed at the world average 0. Table || lists the 
resulting w/ values together with the fraction of the events (//) in each r interval. 

The fcp vertex is determined using lepton tracks from 7/\(f or tl/(25) decays, or prompt 
tracks from Tjc decays. The /tag vertex is determined from well reconstructed tracks not 
assigned to fcp- Tracks that form a Ks are not used. The tracks used for vertexing must 
have at least one three dimensional point with consistent hits in r — (|) plus at least one 
additional z hit. Each vertex position is required to be consistent with the interaction 
point profile smeared in the r-(|) plane by the B meson decay length. We use an iterative 








TABLE 3. The values of sin2(|)i for various 
subsamples (statistical errors only). 


Sample 

sin 2(|)i 

/tag = S'* (? = +l) 

0.84±0.21 

/tag = S'* (^=-1) 

1.11±0.17 

J/^Ks{tC%-) 

0.81 ±0.20 

{cc)Ks except J/MfKs{TC 11^) 

1.00 ±0.40 

J/^Kl 

1.31 ±0.23 


0.85 ±1.45 

All 

0.99±0.14 


procedure: if the quality of the vertex fit is poor, the traek that is the largest eontributor to 
the is removed and the fit is repeated. The typieal vertex-finding effieieney and vertex 
resolution (rms) for zcp (Ztag) are 92 (91)% and 75 (140) ^um, respectively. Note that 
the resolution on the tag side ineludes a large additional eontribution from eharm deeay. 
Onee the tag and CP eigenstate vertiees are reeonstrueted, the proper time is ealeulated 
from Az/y(3. 

After vertexing there are 560 events with q = -\-1 flavor tags and 577 events with 
q = —1. Figure shows the observed Af distributions for the q^f = -|-1 (solid points) 
and q^f — (open points) event samples. In the raw data there is a elear asymmetry 
between the two distributions; this demonstrates visually that CP symmetry is violated. 

To extraet the measured value of the CP violating parameter, we perform an unbinned 
maximum likelihood fit to the time distributions of the tagged and vertexed events. The 
fit takes into aeeount the effeets of baekground, vertex resolutions and ineorreet tagging. 

For modes other than J the pdf expeeted for the signal is 

2’sig(AC^,W;,^/) 

-|Af|/T Q 

= 21^0 {1-2w/)sin2(|)isin(AmrfAf)}, 

where we fix x^o and Am^ at their world average values 0. For the mode a 

more complex fit to Af and the transversity angle[|^ is used to take into aeeount the two 
CP eigenstates whieh eontribute0. 

The pdf used for the baekground distribution is ff’bkg(Af) = + (1 — 

/x) 6( At), where /x is the fraetion of the baekground eomponent with an effeetive lifetime 
Xbkg and 5 is the Dirac delta function. For all fcp modes other than y/\|/.fir/., a study using 
events in the A£' va. Mbc sideband regions shows that the /x component is negligible. For 
these modes we use 2’bkg(At) = 5(At). 

In the ease of 5° — '^Ki, the background is dominated by 5 —> decays where some 
final states are CP eigenstates. We estimate the fractions of the baekground components 
with and without a true Kl eluster by fitting the pg™ distribution to the expeeted shapes 
determined from MC. We also use the MC to determine the fraetion of events with 
definite CP content within eaeh component. The result is a background that is 71% 
from non-CP modes with Xbkg = 'tg. For the CP-mode baekgrounds, we use the signal 
pdf given above with the appropriate values. For whieh is 13% of the 

baekground, we use the '^f = —\ eontent determined from the full sample. The 
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FIGURE 4. The proper decay time distributions for D^*'^ 7t+ decays. The upper curve is the result 
of the fit for the lifetime; the lower curve represents the background contribution. 


remaining backgrounds are — 1 states (10%) including \\fKs, and = +1 states 
(5%) including \\f{2S)Ki, XciKi and \|/7t°. 

The last ingredient needed for the CP fit is the proper-time interval resolution, R{At). 
This resolution function is parameterized by convolving a sum of two Gaussians (a 
main component due to the SVD vertex resolution and charmed meson lifetimes, plus 
a tail component caused by poorly reconstructed tracks) with a function that takes into 
account the cm motion of the B mesons. The relative fraction of the main Gaussian is 
determined to be 0.97 ±0.02 from a study of 5® —> Z)*"p+, D ti^, \\fKs 

and 5+ events. The means (jimam, /^taii) and widths (Omaim cttaii) of the 

Gaussians are calculated event-by-event from the fcp and /tag vertex-fit error matrices 
and the y} values of the fit; typical values are timam = —0.24 ps, /rtaii = 0.18 ps and 
O main = 1-49 ps, Otaii = 3.85 ps. An example of a fit for hadronic non-CP eigenstate 
using the resolution function modes is shown in Fig. ^ The fit agrees well with data out 
to 10 lifetimes on the logarithmic scale of the figure. As a consistency check, we obtain 
lifetimes for the neutral and charged B mesons using the same vertexing procedure that 
is used for the CP fit. The results of the fit agree well with the world average 5® lifetime 
value. 0] 

The pdfs for signal and background are convolved with R{At) to determine the 
likelihood value for each event as a function of sin2(|)i: 

■^i — J At) 

±(1 /sig)±bkg(^^ )Rbkgi^ )}dAt , 

where f^ig is the probability that the event is signal. The most probable sin2(|)i is the 
value that maximizes the likelihood function L = WiLi, where the product is over all 
events. (Note that the signal and background resolution functions are different.) 
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FIGURE 5. (a) The asymmetry obtained from separate fits to each At bin for the full data sample; the 

curve is the result of the global fit. The corresponding plots for the (b) {cc)Ks (^/ = —1), (c) J/\\tKL 
(^f = +1), and (d) control samples are also shown. The curves are the results of the fit applied 
separately to the individual data samples. 


The result of the fit is 

sin2(|)i = 0.99±0.14(stat) ±0.06(syst). 

In Fig. ^a) we show the asymmetries for the eombined data sample that are obtained by 
applying the fit to the events in eaeh At bin separately. The smooth eurve is the result of 
the global unbinned fit. Figures ^b) and (e) show the eorresponding asymmetry displays 
for the {cc)Ks (^/ = —1) and the J/'^Ki = +1) modes separately. The observed 
asymmetries for the different CP states are indeed opposite, as expeeted. The eurves are 
the results of unbinned fits applied separately to the two samples; the resultant sin2(|)i 
values are 0.84±0.17(stat) and 1.31 ±0.23(stat), respeetively. 

The systematie error is dominated by uneertainties in the tails of the vertex distribu¬ 
tions, whieh eontribute 0.04. Other signifieant eontributions eome from uneertainties (a) 





















in wi (0.03); (b) in the parameters of the resolution funetion (0.02); and (e) in the 7 /\|/^l 
baekground fraction (0.02). The errors introduced by uncertainties in Am^ and x^o are 
negligible. 

A number of checks on the measurement were performed. Table lists the results 
obtained by applying the same analysis to various subsamples. All values are statistically 
consistent with each other. The result is unchanged if we use the w/’s determined 
separately for /tag = BP and BP . Fitting to the non-CF eigenstate self-tagged modes 
B^ —> Z)*^p+, J andD*“£+v, where no asymmetry is expected, 

yields 0.05 ± 0.04. The asymmetry distribution for this control sample is shown in 
Fig. Kd). 

A single consistent reconstruction and analysis procedure is used for all data samples. 
We verify that, as expected, the sin2(|)i values for the different run periods are consistent. 
With toy Monte Carlo studies we also verify that the analysis procedure and the errors in 
the analysis are well behaved when sin2(|)i is close to the physical boundary. As a further 
check, we used three independent CP fitting programs and two different algorithms for 
the /tag vertexing and found no discrepancy. 

Finally, we comment on the possibility of direct CP violation. The signal pdf for a 
neutral B meson decaying into a CP eigenstate can be expressed in the more general 
form 




g-|Ar|/T^o l + \XP 

+«,(l-2w,)e,-/m>.sin(A„yA,) 


1-W- 

2 


cos(Am^/Af)]} 


where X is a complex parameter that depends on both 5®^ mixing and on the amplitudes 
for 5® and 5^ decay to a CP eigenstate. The presence of a cosine term (|?i| 7 ^ 1) indicates 
direct CP violation. For the primary analysis, we assumed |?i| = 1 which is the SM 
expectation. In order to test this assumption, we also performed a fit using the above 
expression with /mX (= “sin2(|)i”) and |?i| as free parameters. We obtain |?i| = 1.03 ±0.09 
and ImX = 0.99 ± 0.14 for all CP modes combined, where the errors are statistical only. 
This result confirms the assumption used in our analysis. 


CONCLUSIONS 

In the summer of 2001, Belle (along with BABAR) presented its first significant mea¬ 
surement of the CP violating parameter sin2(|)i. BELLE found 

sin 24)1 =0.99 ±0.14 ±0.06 

with a statistical significance of greater than six standard deviations [PD||. This can be 
compared to the BABAR result of sin2(|)i = 0.59 ±0.14 ± 0.05[[TT||. The two results are 
based on data samples of comparable size (31 million and 32 million BB pairs, respec¬ 
tively). The efficiencies and resolutions of the two experiments are also quite similar. 
However, although the weighted average of the two results agrees well with indirect 
determinations that assume the Standard Model, the two measurements themselves are 





only marginally consistent. Larger data samples and additional more precise measure¬ 
ments will be required to fully reconcile these two results for sin2(|)i. In parallel, a pro¬ 
gram for the measurement of the other angles, (|)2 and (|) 3 , has started. 
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